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ABSTRACT 


An  experimental  study  of  transient  ablation  has  been  carried  out  in 
a  facility  developed  expressLy  for  this  purpose.  The  facility  utiLizes  the 
PIBAL  pebble  bed  convection  heater  as  a  source  of  high  energy  air. 

Hemispherical  models  made  of  severaL  plastic  materials  have  been 
abLated  in  this  facility  and  the  final  configuration  compared  with  the  pre¬ 
dictions  of  a  transient  theory,  which  uses  an  integral  technique.  Good 
agreement  has  been  found  only  for  those  pLastics  which  may  be  classified 
as  of  a  subLiming  type.  For  the  remaining  materiaLu  which  meLt  and  flow 
without  significant  vaporization  the  transient  theory  predicts  vaLues  of  the 
surface  recession  which  are  much  lower  than  those  obtained  experimentally. 
This  discrepancy  many  be  attributable  to  unrealistic  vaLues  of  the  physical 
properties  avaiLabLe  f)r  these  latter  materials. 

The  modeLs  were  exposed  to  both  a  Laminar  and  transttionaL  heat 
transfer  environment.  The  appropriate  heat  transfer  theories  required 
for  the  caLcuLation  are  reviewed  and  discussed  in  detaiL.  Corresponding 
heating  rates  have  been  measured  experimentally  and  are  compared  with 
these  analyses.  Theory  and  experiment  are  shown  to  be  in  good  agreement. 

A  description  of  the  faciLity  and  its  development  is  presented  in 
this  report.  In  addition,  a  survey  of  current  ablation  theories  is  incLuded 
in  the  introduction. 
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I.  INTRODUCTION 


A.  Definition  and  Description  of  the  Ablation  Phenomenon 

The  interaction  between  the  earth's  atmosphere  and  a  solid  body 
traveling  at  hypersonic  speeds  has  been  studied  for  many  years.  Long 
before  the  first  V-2  rocket  was  launched,  astrophysicists  (references  1 
and  2)  were  concerned  with  this  problem  in  connection  with  meteors. 
Recently,  the  aerodynamic  heating  problems  associated  with  the  reentry 
of  ballistic  missiles,  satellites,  and  space  vehicles,  and  with  high  per¬ 
formance  propulsive  systems  have  intensified  the  investigation  of  this  inter¬ 
action.  A  variety  of  heat  protection  systems  have  been  developed  for 
reentering  vehicles.  It  is  interesting  that  one  of  the  most  currently  popular 
methods,  namely  that  of  ablation,  bears  a  strong  resemblance  to  the  meteor 
problem  studied  in  the  past. 

In  its  most  general  form  the  ablation  phenomenon  occurring  during 
reentry  represents  an  exceedingly  complex  problem  in  mathematical 
physics.  Its  solution  requires  a  study  of  inviscid  hypersonic  flow,  the 
heterogeneous  comprensible  boundary  layer,  heat  conduction,  chemical 
reactions,  dissociation  and  ionization,  thermal  stress  analysis,  and  radia¬ 
tion,  as  well  as  the  dynamics  necessary  for  predicting  the  trajectory  of  a 
given  body  with  specific  initial  conditions.  Depending  on  the  point  of  view 
many  of  these  aspects  may  be  treated  independently. 

The  discussion  to  be  presented  here  will  be  somewhat  simplified 
in  form.  Some  aspects  of  the  general  problem  will  be  completely  ignored, 
as  they  properly  belong  to  fields  distinct  from  aerodynamics.  Certain 
additional  information  necessary  for  the  solution  of  the  problem  will  tacitly 
be  assumed  known.*  For  purposes  of  this  discussion,  consider  the 
stagnation  region  of  a  blunt  solid  body  suddenly  exposed  to  a  high  energy 


*  For  example,  the  heat  transfer  rates  to  a  non-ablating  body  of  arbitrary 
shape  are  required  for  certain  aspects  of  the  development.  This  require¬ 
ment  alone  represents  a  major  difficulty  and  continues  to  be  the  subject 
of  intensive  study  at  the  present  time  (see,  for  example,  references  3 
and  4). 

Manuscript  released  for  publication  24  February  1961  as  a  WADD 
Technical  Note. 
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stream  of  air.  Assume  that  the  initial  temperature  within  the  body  ir 
uniform  and  that  the  material  melts  at  an  elevated  temperature.  Until  the 
surface  of  the  body  reaches  the  melting  temperature,  the  heat  which 
diffuses  to  the  surface  will  merely  be  conducted  into  the  interior  of  the 
solid;  during  this  initial  phase  of  heating  the  boundary  layer  and  heat  con¬ 
duction  problems  are  uncoupled  and  may  be  treated  independently.  How¬ 
ever,  as  soon  as  a  phase  change  does  occur  there  arises  a  coupling  of  the 
two  problems,  with  a  resulting  increase  in  complexity. 

As  an  example  consider  the  case  when  the  heating  is  sufficient  to 
cause  vaporization  of  the  solid,  with  or  without  an  intermediate  liquid 
phase.  In  either  case  the  injection  of  gaseous  products  into  the  boundary 
layer  results  in  a  profound  change  in  its  character,  adding  new  constit¬ 
uents  to  the  mixture,  thickening  the  boundary  Layer  and  altering  the  con¬ 
centration,  velocity,  and  enthalpy  profiles.  In  general  the  result  of  this 
mass  transfer  is  a  reduction  of  the  heat  transfer  rates  at  the  interface. 
Obviously,  the  degree  to  which  the  heat  fLuxes  are  reduced  depends  on  the 
rates  at  which  the  gases  are  evoLved  and  vice  versa;  hence  the  coupling 
between  the  two  previously  independent  systems. 

Many  additional  complications  may  arise  when  melting  and  mass 
transfer  occur.  For  example,  if  both  liquid  and  gaseous  phases  coexist, 
the  percentage  of  liquid  which  vaporizes  before  it  is  swept  away  by  aerody¬ 
namic  stresses  must  be  determined.  Of  course,  considerable  simplification 
is  possible  if  the  solid  materiaL  has  a  unique  and  clearly  defined  melting 
and/or  boiling  temperature.  However,  there  are  materials  (the  glasses, 
for  example)  which  are  not  really  solids  but  rather  supercooled  Liquids. 

Such  materials  have  a  very  high  viscosity  even  at  elevated  temperatures, 
so  that  no  distinct  transition  from  solid  to  liquid  phase  is  observed.  The 
temperature  at  the  gas -liquid  interface  becomes  another  unknown,  resulting 
in  an  interdependence  of  the  aerodynamic  and  heat  conduction  problems.  It 
may  be  noted  that  for  such  materials  the  surface  temperatures  may  reach 
values  high  enough  for  a  significant  quantity  of  energy  to  be  radiated  from 
the  surface.  Evidently  such  radiation  is  desirable  as  an  additional  mech¬ 
anism  for  relieving  the  unfavorable  thermal  environment  to  which  the  body 
is  subjected,  provided  that  excessive  heat  transfer  to  the  inner  surface  of 
the  heat  shield  does  not  result  from  the  elevated  surface  temperature. 

The  gaseous  products  injected  into  the  boundary  Layer  may  react 
chemically  with  the  original  constituents  of  the  boundary  layer.  In  general, 
such  a  reaction  can  be  expected  to  be  exothermic,  resulting  in  an  increase 
in  the  heating  potential  to  which  the  surface  is  exposed.  The  chemical 
kinetics  of  the  reaction  must  of  course  be  known  if  its  effect  on  the  overall 
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process  is  to  be  quantitatively  determined. 

The  phase  changes  which  occur  at  the  suriace  of  the  solid  also 
influence  to  a  great  extent  the  manner  in  which  the  heat  conduction  problem 
in  the  interior  of  the  solid  must  be  treated.  The  presence  of  a  liquid  layer 
and  the  evolution  of  gase^  by  evaporation  (or  sublimation)  at  the  liquid-gas 
(or  solid-gas)  interface  introduces  a  nonlinear  character  because  of  the 
presence  of  moving  boundaries.  Such  problems  are  usually  referred  to  as 
‘ ' Stefan -likeM  (reference  5);  exact  analytic  solutions  for  them  are  in 
general  not  available. 

The  presence  of  a  liquid  film  on  the  surface  of  the  solid  results  in 
the  possible  removal  of  material  by  both  the  aerodynamic  forces  as  men¬ 
tioned  previously  and  by  body  forces  which  can  arise  from  rapid  decelera¬ 
tions  of  reentry  bodies.  In  this  connection  the  work  of  Cheng  (reference  6) 
may  be  cited.  It  is  demonstrated  therein  that  rapid  deceleration  can  result 
in  an  unstable  liquid-gas  interface  which  Leads  to  a  direct  Loss  of  mass  as 
Liquid  droplets. 


B.  Review  of  Theoretical  Works  on  AbLation 

The  theoretical  treatment  of  the  exceedingly  complex  ablation 
phenomenon  requires  a  considerable  idealization  of  the  actual  physical 
situation.  In  general,  the  analyses  which  have  been  performed  are  con¬ 
cerned  with  the  simplified  case  of 'a  stagnation  region  on  blunt,  two- 
dimensional  or  axisymmetric  bodies-.  Such  an  approach  permits  the 
assumption  of  laminar  flow  and,  of  course,  gives  the  boundary  layer  equa¬ 
tions  the  "similar"  character  so  essential  to  their  solution. 

For  purposes  of  discussion,  ablation  may  be  categorized  as 
occurring  with  melting  only,  with  sublimation  directly  from  the  solid,  or 
with  a  combination  of  melting  and  evaporation.  In  addition,  one  may  con¬ 
sider  alternately  the  cases  where  the  gaseous  products  either  react  with 
the  constituents  of  the  external  flow  or  are  chemically  inert. 

An  early  work  which  studies  the  case  of  melting  ablation  is  by 
Landau  (reference  7)  who  assumed  total,  instantaneous  removal  of  the 
ablated  material  under  the  condition  of  a  constant  heat  input  at  the  surface. 
Hence  this  study  does  not  consider  the  interaction  of  the  ablation  products 
with  the  external  boundary  layer  and  is  essentially  a  problem  in  heat  con¬ 
duction.  The  theory  predicts  the  melting  rate  and  the  thickness  of  melted 
material  as  functions  of  time. 
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Goodman  (reference  8)  has  considered  the  condition  under  which  the 
assumption  of  instantaneous  removal  of  the  decomposition  products  is  valid 
and  has  shown  inter  alia  that  the  solid  region  must  be  semi-infinite  in  extent 
and  that  its  surface  temperature  be  held  constant  at  its  melting  value. 

Citron  (reference  9)  has  carried  out  a  similar  analysis  for  the  case 
of  a  finite  solid  region  insulated  at  the  unheated  surface  and  has  indicated 
that  the  treatments  of  Landau  and  the  related  work  of  Goodman  cited  above, 
and  of  Sutton  (reference  10)  yield  melting  rates  which  are  too  low. 

Sutton  presented  one  of  the  first  rigorous  analyses  of  melting  abla¬ 
tion  and  its  effect  on  the  external  boundary  layer.  Sutton  has  treated  the 
case  of  Pyrex  glass  subjected  to  hypersonic  flight  conditions  and  integrated 
the  liquid  layer  equations  numerically.  The  importance  of  the  viscosity  of 
the  molten  material  in  blocking  the  heat  conducted  into  the  solid  structure 
is  pointed  out. 

Lees  (reference  11)  considered  the  case  of  melting  ablation  under 
the  assumption  of  constant  viscosity  in  the  liquid  layer  and  obtained  approx¬ 
imate  analytic  solutions  valid  around  body  contours  up  to  the  sonic  point. 

It  was  also  shown  herein  that  the  stagnation  enthalpy  of  the  free -stream  gas 
is  an  essential  parameter  which  must  be  reproduced  in  experimental  work 
if  similitude  is  to  be  achieved. 

Bethe  and  Adams  (reference  12)  treat  the  ablation  of  a  glassy 
material  (i.  e.  ,  one  with  an  exponential  variation  of  viscosity  with  temper¬ 
ature)  at  a  stagnation  region,  and  consider  the  effect  of  partial  vaporization 
of  the  molten  film.  There  are  obtained  approximate  solutions  which  deter¬ 
mine  the  energy  absorbed  per  unit  mass  of  ablated  material.  In  addition, 
the  dependence  of  this  and  other  ablation  parameters  on  flight  conditions  is 
established. 

The  case  of  subliming  ablation,  which  occurs  with  some  plastic 
materials,  is  particularly  well  suited  to  treatment  by  a  simplified  engi¬ 
neering  approach.  Typically,  it  is  assumed  that  the  material  has  a  unique 
sublimation  or  ablation  temperature  which  is  constant,  thai  the  phenomenon 
is  qua3i-steady  (i.e.  ,  ablation  proceeds  at  a  constant  rate),  and  that  the 
thermal  conductivity  is  small,  so  that  the  thermal  penetration  into  the  solid 
is  negligibly  small.  The  work  of  Roberts  (reference  13)  and  Georgiev, 
Hidalgo,  and  Adams  (reference  14)  can  be  cited  in  this  connection.  The 
concept  of  the  effective  heat  of  ablation  (hg££  )  is  utilized  therein  and  esti¬ 
mates  of  its  magnitude  and  its  dependence  on  environmental  conditions 
obtained  by  the  analysis.  Some  experimental  results  are  presented  in 
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reference  14.  Additional  data  may  be  found  in  references  15  and  16. 

Much  of  the  analysis  of  ablation  with  melting  and  evaporation  has 
been  carried  out  with  reference  to  certain  specialized  types  of  materials. 

In  particular,  these  are  the  gLas s-plastic  combinations  and  the  high- 
viscosity  glasses  such  as  quartz  and  Pyrex.  The  former  materials,  which 
may  also  be  referred  to  as  reinforced  plastics,  are  produced  by  molding 
a  mixture  of  a  thermosetting  plastic  and  refractory  fibers  (such  as 
fiberglas),  the  plastic  serving  as  the  matrix  for  the  reinforcing  fibers. 

The  concept  of  a  gasification  ratio  F  introduced  by  Scala  (references  17  and 
18)  is  useful  in  this  case  and  represents  an  additional  unknown  of  the 
problem.* 

Typical  of  the  theoretical  and  experimental  investigations  for  a 
composite  materiaL  is  the  work  ofGeorgiev,  Hidalgo,  and  Adams  (refer¬ 
ence  19)  who  consider  the  case  of  glass-reinforced  phenolic  plastic.  It  was 
found  that  the  gLass  flowed  without  evaporation  while  the  plastic  matrix 
sublimed.  The  resulting  interface  temperature  was  shown  to  depend  on  the 
liquid  flow,  i.  e. ,  to  be  determined  by  the  solution  o£  the  liquid  layer  equa¬ 
tions.  Additional  work  by  Scala  may  be  cited  (reference  20)  wherein  a 
dependence  of  heff  and  T  on  stagnation  pressure  as  well  as  on  stagnation 
enthalpy  is  deduced. 

Representative  anaLyses  for  the  glassy  materials  may  be  found  in 
references  21  and  22.  In  the  former  a  momentum  integral  technique  is 
utilized  for  the  determination  of  the  steady-state  ablation  characteristics  of 
the  liquid  iayer.  These  equations,  with  assumed  forms  for  the  temperature 
and  velocity  distribution  in  the  liquid  regime,  reduce  to  aLgebraic  equations 
at  the  stagnation  region.  These  are  then  matched  with  exact  solutions  for 
the  dissociated  gas  phase  (reference  23)  to  determine  the  various  unknowns 
of  the  problem.  Details  of  this  procedure  are  given  in  reference  24.  In 
reference  22,  the  identical  probLexn  is  considered  with  the  liquid  assumed 
to  be  incompressible  and  with  a  temperature  dependent  viscosity.  The 
governing  equations  are  dedeuced  and  temperature,  shear  stresses,  mass 
and  energy  transfer  matched  at  the  gas-liquid,  liquid-solid  interfaces.  A 
numerical  solution  for  Pyrex  gLass  is  presented. 

An  experimental  and  theoretical  study  in  which  the  capabilities  of 
the  various  types  of  ablation  are  compared,  may  be  found  in  references  18 
and  21.  It  is  shown  therein  that  the  selection  of  the  "best"  materiaL 


The  gasification  ratio  T  represents  the  fractional  part  of  the  solid  which 
gasifies. 
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depends  critically  on  fLight  conditions,  that  is,  on  the  type  o£  application 
intended.  This  latter  aspect  will  be  discussed  briefly  in  the  next  section. 

Mac  C.  Adams  (reference  25),  in  an  extensive  ar.d  thorough  survey 
article,  presents  correlation  formulas  for  the  effective  heat  of  ablation 
for  quartz  and  for  a  glass-plastic  combination.  These  relations  are 
based  on  the  work  performed  in  some  of  the  aforementioned  references 
(11,  19,  and  ?.6)  and  yield  simple  engineering  formulas  which  relate  heff 
to  free-stream  stagnation  pressure  and  enthalpy.  Estimates  for  the  case 
of  turbulent  heating  are  aLso  indicated. 

Ablation  with  chemical  reaction  has  been  treated  by  several 
investigators  (references  18,  27,  and  28).  In  particular,  the  work  of  Knuth 
(reference  27)  and  Denison  and  Dooley  (reference  28)  will  be  outlined  here 
as  being  representative.  Any  analysis  of  this  problem  requires  rather 
detailed  information  concerning  the  chemical  kinetics  of  the  ablating 
material.  In  general,  such  information  is  not  available  for  technologically 
interesting  reactions  so  that  at  present  quantitative  results  can  only  be 
obtained  for  a  few  materials.  The  references  cited  above  consider  the 
ablation  (with  combustion)  of  Teflon  and  graphite. 

Knuth  treated  the  steady  Couette  flow  with  diffusion  of  a  reacting 
gas  (the  Teflon  monomer  CaF4).  Under  the  assumption  of  infinitely  fast 
(and  infinitely  slow)  reaction  rates  the  pertinent  equations  are  deduced, 
and  velocity,  enthalpy,  and  partial  pressure  profiles  calculated  therefrom. 
An  important  result  is  that  the  ablation  rate  of  Teflon  when  exposed  to  a 
hypersonic  environment  is  influenced  insignificantly  by  combustion  of 
the  monomer. 

In  reference  28,  a  detailed  aerothermochemical  analysis  of  chem¬ 
ically  active  surfaces  in  a  high  speed  airstream  of  an  oxidizing  gas  mix¬ 
ture  is  presented  and  applied  to  the  particular  case  of  graphite.  In 
agreement  with  several  other  works  (references  27,  29,  and  30),  it  is 
shown  that  the  exothermic  reaction  does  not  greatly  affect  the  mass  and 
heat  transfer  at  the  solid  surface  when  the  heating  rate  is  sufficiently 
high  to  cause  the  reaction  to  occur  within  the  boundary  layer. 

C.  Engineering  Significance  and  Applications  of  Ablation 

The  effectiveness  of  ablation  in  extreme  thermochemical  environ¬ 
ments  is  most  evident  when  it  is  compared  with  other  methods  of  pro¬ 
tection,  such  as  the  heat  sink,  radiators,  and  transpiration  cooling.  A 
comparison  of  these  on  a  weight  basis  indicates  that  ablation  is  one  form 
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or  another  (see  previous  remarks  about  various  types  of  ablation)  is  the 
most  practical  and  effective  for  virtually  all  of  the  aforementioned  environ¬ 
mental  conditions  (see,  for  example,  reference  17).  Another  advantage  of 
the  ablation  is  its  passivity  and  the  self- regulation  which  it  exhibits;  that  is, 
as  the  heat  transfer  rate  to  the  ablating  surface  increases,  the  rate  of 
ablation  increases  so  as  to  absorb  the  increase  and  to  reduce  the  heat 
transfer  rate.  The  practicability  of  the  abLation  process  in  missile  appli¬ 
cations  has  already  been  established. 


D.  Summary 

The  references  cited  above  are  representative  of  the  effort  devoted 
recently  to  the  ablation  phenomenon.  Although  a  great  part  of  this  work 
utilizes  a  highly  theoretical  approach  which  yields  a  broader  understanding 
of  the  underlying  principles  governing  the  process,  it  is  interesting  to  note 
that  the  major  contribution  to  the  solution  of  the  reentry  problem  has  been 
due  to  research  involving  an  elementary  engineering  point  of  view  in  com¬ 
bination  with  experiment.  The  major  factor  which  distinguishes  the  two 
approaches  is  that  in  the  latter  case  steady-state  ablation  is  considered. 

It  is  evident,  however,  that  in  order  to  make  final  design  estimates  and  to 
check  the  validity  of  these  approximate  calculations,  the  more  complicated 
transient  ablation  analysis  must  be  carried  out.  Consequently,  experimental 
techniques  for  the  study  of  the  transient  case  must  also  be  developed. 

For  materials  of  practical  interest  exposed  to  the  high  heat  fluxes 
corresponding  to  typical  reentry  conditions,  the  accurate  measurement  of 
the  transient  behavior  is  exceedingly  difficult,  due  to  the  relatively  rapid 
rate  with  which  the  abLation  becomes  steady.  One  way  by  which  this  diffi¬ 
culty  can  be  overcome  is  by  exposing  the  modeLs  to  air  streams  of  relatively 
low  stagnation  enthalpies.  In  this  way  the  transient  portion  of  the  process 
is  extended  in  time  so  that  accurate  measurements  may  be  obtained.  If  in 
addition  the  models  are  large,  further  precision  results. 

The  facility  and  the  experimental  techniques  which  are  described 
in  this  report  were  developed  in  keeping  with  this  point  of  view.  The  wind 
tunnel  used  to  ablate  the  models  which  were  tested  produces  an  airstream 
capable  of  ablating  the  specimens  with  a  significant  transient  period;  this 
permits  a  photographic  resolution  with  time  of  the  unsteady  surface  reces¬ 
sion.  The  size  of  the  models  which  were  designed  enhanced  the  precision 
of  these  measurements  and  also  permitted  the  installation  of  thermocouples 
below  the  surface,  so  that  additional  data  of  interest  could  be  obtained.  In 
the  next  section  the  theoretical  consideration  required  for  the  reduction  of 
the  data  is  presented.  Then  follows  a  description  of  the  development  of  the 
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facility  and  the  design  of  the  models  and  test  procedure  utilized.  Finally, 
the  data  whichhave  been  obtained  are  discussed  and  compared  with  the 
pertinent  theories. 


II.  THEORETICAL  ANALYSES 


A.  Transient  Ablation 

The  transient  sublimation  of  a  semi-infinite,  one -dimensional  solid 
with  constant  physical  properties,  exposed  to  a  hypersonic  airstream  with 
constant  fLow  properties,  is  considered.  In  particular,  it  is  assumed  that 
the  solid  sublimes  at  a  definite  and  unique  temperature  hereinafter  referred 
to  as  the  sublimation  temperature  and  denoted  by  T8.  The  analysis  is 
generally  specialized  to  the  conditions  under  which  the  experimental  portion 
of  this  investigation  has  been  conducted. 

The  phenomenon  of  transient  sublimation  can  be  separated  into 
three  distinct  time  intervals,  each  involving  a  different  mathematical  des¬ 
cription  or  analysis.  These  time  intervals  can  be  defined  as 


0  <  t  <  tB 


t  <t 


ss 


t  <  t  <  GO 

SS  —  — 


where  t8  is  the  time  required  for  the  surface  of  the  solid  material  to  reach 
the  sublimation  temperature  T8,  and  t88  represents  the  time  required  to 
achieve  a  condition  of  steady  state  sublimation.  Mathematically,  this 
latter  time  is  indistinguishable  from  t — >oo,  but  from  a  physical  point  of 
view  this  distinction  is  warranted,  as  will  be  seen  below.  The  time  t=0 
corresponds  to  the  instant  of  application  of  the  heat  flux. 

Initially  the  solid  is  assumed  to  be  at  a  uniform  temperature 
Ti(Ti  <Ta).  At  t=0  it  is  instantaneously  exposed  to  a  hypersonic  airstream 
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with  steady  flow  properties.  UntiL  the  surface  of  the  solid  reaches  the 
temperature  T»,  the  solid  can  be  considered  to  be  an  idealized  heat  sink 
with  the  temperature -time  history  described  by  the  unsteady,  one¬ 
dimensional  heat  conduction  equation  (see  Fig.  1  for  coordinate  system) 


..  8  T(x,t) 

x  ■  ■  1  » - 

8x 


8T 

dt 


0  <  t  <  t8 


(1) 


with  initial  condition 


T(x,0)  =  Tt 


(2) 


and  boundary  conditions 


Llm 


oo 


T  (x,t)  = 


(3) 


k  =  h[Too-T<o.t)] 


(4) 


Here.  H  and  k  represent  the  thermal  diffusivity  and  conductivity  of  the 
solid,  respectively,  and  h  tho  heat  transfer  or  film  coefficient  defined  by 
the  relation 


h  = 


TQO-T(0,t) 


(5) 


and  Too  ls 


where  qQ  is  the  net  heat  flux  per  unit  area  to  the  surface  x=0, 
the  stagnation  temperature  of  the  external  flow.  The  film  coefficient  h  is 
assumed  to  be  constant  in  this  analysis.  The  solution  to  the  system 
(l)-(4)  can  be  found  in  the  literature  (see,  for  example,  reference  31)  and 
is  given  by 


T (x,t)  =  \  +  (Too-Ti)  [erfc  X-eNu(Nu+2X)  erfc  (Nu+X)] 


(6) 
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where 


X  = 


2  */kT 


(7) 


the  notation  being  that  of  reference  32. 

The  quantity  t8  can  now  be  evaluated  from  the  relation 

Nu* 

T  (0,t8)  ■  Ts  =  Ti  +  (T^-Tj)  (l-e1NU  erfc  Nu)  (8) 


obtained  by  setting  X=0  in  Eq.  (6).  The  temperature  distribution  at  time 
t=t8,  namely  T(xtts)f  can  then  be  obtained,  again  utilizing  Eq.  (6). 

When  the  surface  of  the  solid  is  at  the  sublimation  temperature  the 
solid  will  vaporize,  and  the  vapor  will  be  injected  into  the  boundary  Layer, 
where  it  is  removed  by  convection.  The  surface  of  the  solid  will  recede. 
Let  x~s( t)  denote  the  distance  of  this  receding  surface  from  its  initial 
position  x=0.  Thus,  x=s(t)  represents  a  moving  interface  between  the  solid 
material  and  the  external  gas  boundary  Layer.  If  q(t)  is  the  net  heat  flux  to 
this  surface,  L  the  heat  of  subLlmation  of  the  solid  and  p  its  density,  an 
energy  balance  at  the  interface  requires  that 


k 

8x 


ds 

q(t)  -  PL  — 
dt 


(9) 


That  is,  the  net  heat  flux  to  the  solid  is  partLy  utilized  to  sublime  some  of 
the  soLid  while  the  remainder  is  conducted  into  the  interior. 

For  the  process  described  above  an  unsteady  heat  conduction 
problem  can  be  formulated  as  follows:  for  s  <  x  <  oo  and  tg  <  t,  the  heat 
conduction  Eq.  (1)  is  valid  and  the  following  boundary  conditions  apply: 


lim  8T(x,t)  =  Q 
x — >oo  8x 


(10) 
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(11) 


T(e,t)  =  Tb 


while  the  initial  spatial  temperature  distribution  is  given  by  the  constant 
h  analysis  presented  previously,  namely: 


T(x,ts) 


(12) 


The  solution  to  the  problem  defined  by  the  system  of  Eqs.  (1),  (10), 
(11),  and  (12)  presents  considerable  difficulty  due  to  the  existence  of  the 
moving  boundary  s(t).  A  method  analogous  to  the  Karman  and  Polhausen 
momentum  integral  (reference  33)  can  be  used  to  satisfy  the  heat  conduction 
equation  on  the  average  and  to  obtain  an  approximate  but  simple  solution. 
The  "heat-balancc  integral"  thus  obtained  was  introduced  by  Goodman 
(reference  34)  and  applied  to  various  heat  conduction  problems  involving  a 
change  of  phase.  Such  a  technique  does  not  yield  exact  results,  but  simple 
analytic  solutions  can  be  obtained  which  appear  to  have  sufficient  accuracy 
for  most  engineering  purposes. 

The  analysis  presented  below  differs  from  that  in  reference  34  in 
two  important  respects.  First,  in  the  choice  of  a  temperature  profile 
within  the  solid,  and  second  in  the  assumed  form  for  the  heat  flux  rate.  In 
particular,  it  may  be  noted  that  in  the  aforementioned  reference  T(x,t)  is 
taken  to  be  parabolic  in  form,  while  q(t)  is  assumed  to  be  constant.  In 
what  follows  these  are  replaced  by  Eqs.  (13)  and  (14)  which  appear  to  be 
somewhat  more  realistic.  The  justification  for  these  particular  choices  is 
discussed  in  Appendices  I  and  n. 

Let  the  temperature  distribution  within  the  unablated  solid  be  given 
by 


T(x,t)  =  T.  +  (Ts-T.)  exp[- 


x-s(t). 

«(t)  J 


(13) 


In  addition  let 


q(f)  =  qo  ■  n  ar 


(14) 


where  qQ  and  II  will  depend  only  on  the  physical  properties  of  the  solid  and 
on  free-stream  conditions;  their  significance  as  well  as  that  of  the  quantity 
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(t)  Introduced  In  Eq.  (13)  are  discussed  in  the  appendices. 

Eq,  (1)  is  now  integrated  over  the  remaining  soLid,  i.  e. : 


from  which 


x[. 


8T(oott)  9T(fl, 


8x 


8x 


T(x,t)dx  ♦  T(s,t)  £  US) 

dt 


From  the  boundary  conditions  (10)  and  (11)  and  the  assumed  profile 
(13),  Eq.  (15)  becomes 


H  d6  ds 
6  dt  dt 


(16) 


while  (9)  becomes  [utilizing  Eq.  (14)] 


^=qo-(PL+n)  £ 


V  5  (Tg-Tj) 


(17) 


Eliminating  ds/dt  between  Eqs.  (16)  and  (17)  there  is  obtained  a  differential 
equation  for  the  quantity  6(t)  as  follows: 


d8  , 


dt 


kV 

PL+n 


1  + 


qo 

PL+n 


=  o 


(18) 


Eqs.  (17)  and  (18)  represent  two  simultaneous,  ordinary  differential 
equations  for  8  and  s.  The  initial  conditions  for  this  system  are 
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for  t=t8  , 


.  da 

i=0  ,  —  =  0 

dt 


5=S 


Introducing  now  the  transformations 


T  = 


„  *lo 
*n  =  -si.  s 

kV 


,  kV  1 


k(pl+ii) 


the  systems  (17)  and  (18)  become 


^  +  (i+x)  ca-  xr  =  0 


and  the  initial  conditions  transform  as 


T  =  0;  r\  =  0,  £  =  1 


The  solution  to  this  system  is  quite  easily  obtained.  It  is 


t  =  i±l  ln[ - i - J  .  I  (1  .!) 

xa  (1+X)£-XJ  X  i 
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for  X^O.  Note  that  for  Large  times  (i.e.  ,  T 


oo) 


£ 


X 

1  +  X 


T| 


X 

1  +  X 


(26) 


These  conditions  represent  the  steady-state  solution  referred  to 
previously;  that  is,  for  large  times  the  ablation  velocity  is  constant,  as  is 
the  temperature  distribution  within  the  remaining  solid;  this  may  be 
verified  by  inspection  of  the  temperature  profiLe  of  Eq.  (13)  under  the  steady 
state  requirements  outlined  in  Eq.  (26).  The  profile  takes  the  form 


T(x,t)-T. 

—  "T  =  exp[-A(x-Bt+C)]  (27) 

Ts-1i 


With  respect  to  a  coordinate  system  fixed  with  the  steadily  receding 
surface  (i.  e. ,  with  velocity  B)  the  profiles  have  the  form 


T(x',t)-T. 

~  r  exp( -Ax  ) 


T.-Ti 


x  =  x 


Bt  +  C 


(28) 


and  so  the  profiles  are  time  independent.  The  constants  introduced  above 
are  defined  by 
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(29) 


A  = 


B  = 


C  = 


B 

K 


*  v* 

1  +  X  kV 


(the  steady -state  ablation  velocity) 


18  ^  X« 


Eqs.  (24)  and  (25)  are  plotted  in  Fig.  3  for  various  values  of  the 
parameter  X.  Additional  plots  of  the  physical  variables  are  presented  in 
Figs.  4  to  7,  with  Teflon  taken  as  the  subliming  material. 

Consider  now  modification  of  the  above  analysis  to  the  case  of  a 
melting  material.  In  this  case,  consistent  with  the  assumption  of  instan¬ 
taneous  removal  of  the  melted  material,  Eq.  (14)  is  replaced  by  q(t)=q0 
which  is  mathematically  equivalent  to  setting  11=0  in  Eq.  (14)  and  in  the 
definition  of  X,  Eq.  (20).  The  remaining  formulation  is  identical.  Evidently 
Tfi  and  L  are  now  interpreted  as  the  melting  temperature  and  heat  of  fusion, 
respectively.  The  corresponding  time  required  for  the  surface  to  reach 
the  melting  temperature  will  be  denoted  by  tm. 

It  is  of  interest  to  compare  the  results  for  the  melting  material  wjth 
the  analysis  of  Landau  (reference  7).  With  the  exception  of  the  boundary 
condition  which  is  assumed  to  apply  before  melting  starts,  both  analyses 
treat  an  identical  problem.*  Comparison  of  the  two  theories  indicates  that 
the  expressions  derived  for  the  steady  state  ablation  velocity  [Eq.  (29)]  and 
the  steady  state  temperature  profile  [Eq.  (28)]  are  identical  with  the  results 
of  reference  7.  The  amount  of  material  which  is  removed,  as  predicted  by 
the  two  theories,  is  shown  in  Fig.  8  in  terms  of  the  nondimensional  variables 
of  reference  7.  We  note  that  the  agreement  is  good  for  physical  times 
corresponding  to  test  conditions.  The  large  deviation  which  occurs  at  very 


*In  reference  7  it  is  assumed  that  for  all  t>0,  q(t)  is  constant  at  the  surface, 
while  in  the  analysis  here  the  boundary  condition  at  the  surface  is  taken  as 
h=q(t)/[Too-T(0,t)]=constant.  Evidently  for  t  >tm  these  conditions  are 
equivalent,  since  it  is  assumed  that  the  surface  remains  at  this  temper¬ 
ature  thereafter.  The  two  problems  will  therefore  differ  solely  in  the 
manner  in  which  tm  is  evaluated.  For  the  conditions  under  which  the  exper¬ 
imental  portion  of  this  report  was  performed  this  difference  has  a  negligible 
effect  on  numerical  results. 
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FIG.  3  LOCATION  OF  VAPORIZATION  LINE  AND  THERMAL 
PENETRATION  (NON-DIMENSIONAL) 
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FIG.  4  LOCATION  OF  VAPORIZATION  LINE  VS  TIME 

MODEL  T-IO .  s-0 
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FIG.  5  TEMPERATURE  HISTORY  AT  SPECIFIC  DEPTH 
MODEL  TIO  -  TRANSIENT  THEORY 


FIG. 


FIG.  8  COMPARISON  OF  PRESENT  THEORY  WITH 
RESULTS  OF  REFERENCE  7 
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short  times  is  due  to  the  previously  mentioned  difference  in  the  vaLues  of 


t 


m* 


As  has  been  previously  noted  in  the  introduction,  the  assumption  of 
instantaneous  removal  of  the  melted  material  effectively  neglects  the  inter¬ 
action  between  this  Liquid  and  the  external  fLow.  More  detailed  analysis 
reveaLs  that  this  interaction  tends  to  reduce  the  interface  heat  fLux  resulting 
in  a  decrease  in  the  rate  at  which  materiaL  is  removed  (cf.  reference  10). 
Thus  the  more  elementary  theory  considered  here  will  predict  higher 
values  for  the  ablation  rates,  and  in  this  sense  it  is  conservative.  This 
point  wiLL  be  considered  again  in  connection  with  the  experimental  results 
which  were  obtained. 

FinaLly,  it  is  interesting  to  note  that  one  may  deduce  the  thermal 
diffusivity  x  of  the  pLastic  materials  experimentally  by  utilizing  the  analysis 
presented  above.  The  technique  is  due  to  Sutton  (reference  35)  and  will  be 
described  below. 

Consider  the  expression  which  has  been  deduced  for  the  steady 
temperature  distribution,  Eq.  (28): 


T(x  ,t)  -  T£  =  (Ts-Ti)  exp(-^x  ) 


where  B  is  the  steady-state  ablation  velocity.  Here  x"  is  the  distance  from 
the  subliming  surface  and  is  given  in  the  steady  condition  by 


x'=x-Bt+C 


If  we  consider  the  temperature  history  at  a  point  initially  ^  distant  from 
the  surface  we  have  for  the  temperature  history  at  this  Location 

T (t)  -  T.  =  (T,-T.)  exp[-  §•  (xQ-Bt+C)] 


Taking  the  logarithm  of  this  equality,  we  obtain 

ln[T(t)-Ti]  =  D  +  Sit 


(30) 
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where 


D  ■  In  V  -  5.  (Xq  +  C) 


Eq.  (30)  will  plot  as  a  straight  line  with  slope  equal  to  (see  Fig.  6). 
Hence  the  semi-log  plot  of  the  measured  temperature  history  combined 
with  measurements  of  the  steady-state  abLation  velocity,  yields  the  experi¬ 
mentally  the  thermal  diffusivity  of  materials. 


B  .  Steady  Ablation 

The  analyses  to  be  presented  here  is  due  to  S.  Georgiev,  et  al 
(reference  14)  and  is  included  in  order  to  emphasize  some  additional 
parameters  which  are  pertinent  to  the  study  of  ablation  from  an  engineering 
point  of  view. 

It  has  been  shown  in  the  previous  analysis  that  when  a  steady  state 
condition  of  ablation  has  been  achieved  the  temperature  profiles  within  the 
solid  remain  invariant  with  respect  to  the  receding  surface.  In  this  case 
the  energy  within  the  solid  remains  constant.  During  this  period,  there¬ 
fore,  the  heat  conduction  term  may  be  eliminated  from  the  heat  balance  at 
the  ablating  surface  [Eq.  (9)]  with  the  result 


q  -  qr  =  m(cpV+L)  (31) 

Here  m  is  the  steady-state  mass  loss  rate  (m=pB),  a  constant,  and  the 
remaining  symbols  retain  their  previous  meaning.  Eq.  (31)  states  that  the 
aerodynamic  heat  transfer  to  the  solid  is  absorbed  by  the  heat  capacity 
plus  the  heat  of  sublimation  of  the  ablating  material.  It  will  be  noted  that 
an  additional  term  qr  has  been  included  to  account  for  heat  radiation  from 
the  surface. 

The  dependence  of  the  aerodynamic  heat  transfer  rate  on  the  injec¬ 
tion  rate  introduced  previously  [see  Eqs.  (14)  and  (A-3)]  is  again  valid  and 
yields,  together  with  Eq.  (31): 


x 

cpV+L+  ^‘(hoo-h,,) 

q<3  =  mt - ; - 7 - ] 

r  o 


(32) 
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If  now  an  effective  heat  of  ablation  ie  defined  by  ^^sq^m, 


leff 


1 

cpV  +  L+  jM4(h00-h8) 


(33) 


Eq.  (33)  can  be  modified  to  account  roughly  for  the  effect  of  any 
chemical  reaction  of  the  injected  vapor  with  the  external  flow.  In  refer¬ 
ence  36,  with  the  assumptions  Lewis  number  equal  to  unity,  and  total 
reaction  of  all  the  oxygen  diffusing  to  the  surface,  an  expression  for  the 
increase  in  interface  heat  transfer  is  deduced: 


Aq  = 


0.21  h. 


hCO-hB 


q 


(34) 


Here  hc  is  the  heat  released  per  unit  mass  of  oxygen  entering  into  the 
reaction,  and  the  numerical  factor  is  the  mass  fraction  of  Oa  external  to 
the  boundary  layer.  Substitution  of  Eq.  (34)  in  Eq.  (33)  yields 


h 


eff 


1+1. 


1 

|  M*(h 


co 


-h.)]  t1- 


.-1 


q  (i+q  > 

o  c 


(35) 


where 


<lc  5 


0.21  h 


h  -h 
oo 


Z 

s 


It  is  noted  that  Eq.  (35)  predicts  no  dependence  of  the  effective 
enthalpy  on  the  heat  flux  rate  in  the  absence  of  radiation.  In  addition,  for 
sufficiently  high  values  of  the  boundary  layer  enthalpy  difference,  the 
contribution  of  combustion  is  also  negligible,  as  has  been  indicated  prev¬ 
iously  (cf,  reference  27). 

If  the  effect  of  radiation  and  chemical  reaction  are  neglected,  there 
results  a  simple  linear  relation  between  heff  and  the  "driving  enthalpy" 
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(36) 


heff  -  V  +  L  +  !  “  (hoo-h.' 


For  a  given  material,  Cp,  V,  L,  and  M  are  known  quantities  so  that  a 
simple  calculation  immediately  yields  the  heat  absorbing  capability  of  a 
given  material  under  specified  flow  conditions.  Conversely,  a  rather 
simple  experimental  procedure  (i.  e.  ,  measurement  of  qo  and  m  for 
various  values  of  h^-hg)  permits  the  determination  of  the  intercept 
(cpV+L)  and  the  slope  (2/3  of  Eq.  (36),  and  thus  yields  important 

information  concerning  the  properties  of  the  material. 

The  results  described  above  have  been  verified  experimentally. 


C.  Aerodynamic  Heating 

Under  the  conditions  to  which  the  ablating  models  were  exposed  in 
this  investigation  both  laminar  and  transitionaL  heating  may  be  expected. 
Hence,  for  purposes  of  comparison  a  theoretical  analysis  for  both  types 
of  heat  transfer  is  required. 

The  laminar  theory  utilized  is  due  to  Lees  (reference  3)  with 
modification  for  the  more  accurate  stagnation  point  theory  of  Fay,  Riddell, 
and  Kemp  (reference  37).  In  terms  of  the  nondimens ional  parameters 
defined  in  the  List  of  Symbols  this  analysis  gives 


where 


pw  ^w 
P  H 
*e  ‘e 


0.1 

) 


F(s) 


F(s) 


P  IX  u  r 


_ a_! 

P  U  u  r  ds] 


(37) 


(38) 


and  the  Prandtl  number  is  taken  to  be  0.71,  If  an  isentropic  relation 
between  velocity  and  pressure  is  assumed,  then,  for  a  given  pressure  dis¬ 
tribution  and  body  geometry  F(T)  may  be  evaluated,  using  an  average 
isentropic  exponent  together  with  isentropic  flow  tables,  such  as  those 
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given  in  reference  38.  For  specified  external  flow  conditions  the  heat 
transfer  rate  q^,  may  be  computed  from  the  definition 


NNu  =  qoCpt  Ro/kte<Vhw> 


(39) 


For  comparison  with  the  experimental  heat  transfer  data  which  were 
obtained,  tjQ  is  evaluated  from  Eq.  (39)  for  the  wall  temperature  equal  to 

the  initial  surface  temperature  of  the  model;  that  is,  qJ  .  On  the 

IT  w=”i 

other  hand,  the  quantity  which  is  required  for  the  calculation  of  the  perti¬ 
nent  ablation  parameters  is  the  value  of  the  aerodynamic  heat  transfer 
when  the  walL  temperature  is  equal  to  the  ablation  temperature,  namely 
This  latter  quantity  is  obtained  by  recalling  that  the  heat 


Hr 


=T 

w  s 

transfer  coefficient  h  is  constant  in  the  test  section;  then 


s 


Tw=Ti 


=  MT^-T.)  =  (T^-T.) 

Aoo  Ai 


(40) 


It  should  be  noted  that  the  value  of  h  varies  locally  from  point  to  point  on 
the  body  surface  but  is  a  constant  at  any  given  station  T . 

For  the  calculation  of  transitional  heat  transfer  a  numerical 
integration  of  the  momentum  equation  for  the  axi symmetric  boundary  Layer 
was  required.  This  equation  may  be  written  in  the  form 


dNo 

^  =  Nr, 


42-h, 


lln 


LL  \ 


(41) 


where 


e 
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fCO 

b  -  J  (Pu/Peue)(l-u/ue)dn 


the  momentum  thickness. 


n  =  coordinate  normal  to  the  wall,  and 
Cf  =  2  Tw/ Peuea  ,  the  local  skin  friction  coefficient. 


For  the  integration  of  Eq.  (41)  the  relation  between  the  skin  friction 
coefficient  and  must  be  known.  In  accordance  with  the  suggestion  of 

Persh  (reference  39)  and  including  the  modification  of  reference  40,  the 
skin  friction  law  used  here  was 


H 

2 


_  i 

=  iT[0.013N^'5.KNR4-*] 


(42) 


where  K  is  selected  such  that  Cf/2  is  continuous  at  transition.  In  refer¬ 
ence  39,  on  the  basis  of  experimental  data  from  flat  plate  transition,  it  ie 
concluded  that  a=2.  The  results  here  suggest  that  in  a  favorable  pressure 
gradient  the  transition  zone  may  be  more  extended;  accordingly,  a  value 
of  a=l  was  assumed.  The  initial  value  of  NR  at  the  start  of  transition  is 
obtained  from  the  relation 


Nro 


0.66[  f  p  u  r  ds  ]z 

-/ft 


(43) 


Nt 


r  H 


This  expression  corresponds  to  the  laminar  theory  of  Lees  cited  prev¬ 
iously.  Eq.  (41)  may  now  be  integrated  from  the  transition  point  (deter¬ 
mined  experimentally)  and  yields  the  skin  friction  and  momentum  thickness 
distributions.  For  the  calculation  of  the  corresponding  heat  transfer  dis¬ 
tribution,  Reynolds  analogy  (valid  for  the  "cold  wall"  case;  cf.  reference 
40),  modified  to  account  empirically  for  ncnunity  Prandtl  number  effects, 
is  used.  There  results 
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(44) 


-  i 

%  =  ff  *  (haw-hw>  PeUe(cf/2) 

where  CX  =  Prandtl  number  and  haw  is  the  adiabatic  wall  enthalpy  defined 
by  the  equation 


(45) 


In  terms  of  the  Nusselt  number  Nnu ,  Eq.  (44)  may  be  written 


NNu  =  °*89(cf/2)GNRP  u  (46> 


where  G  *  (haw-hw)/(hg  -bw)  and  the  Prandtl  number  is  again  taken  as 
0.71.  Integration  of  Eq.^(41)  was  carried  out  numerically  on  the  PIBAL 
Bendix  G-15  General  Purpose  digital  computer  by  the  Runge-Kutta  method. 

For  fully  developed  turbulent  heat  transfer  the  prediction  of  the 
flat  plate  reference  enthalpy  method  (cf.  Method  I  of  reference  41)  has 
been  used. 


D.  General  Remarks 

It  is  evident  from  the  foregoing  that  the  fundamental  measurements 
to  be  obtained  for  the  study  of  steady  subliming  ablation  are  the  quantities 
q0  and  m.  For  the  transient  case  a  time  resolved  history  of  the  surface 
recession  s(t)  would  be  required.  It  wouLd  also  be  of  interest  to  obtain  the 
temperature -time  history  T(Xo,t)  at  some  point  below  the  surface  of  the 
ablating  material  for  comparison  with  the  transient  theory  which  has  been 
developed.  These  quantities  should,  of  course,  be  measured  over  a  range 
of  stagnation  conditions,  so  that  their  variation  with  hgQ-hg  and  qo  may  be 
established.  Certain  additional  information  such  as  stagnation  temperature 
and  pressure  as  weLL  as  the  existing  pressure  distribution  is  also  required. 

In  the  following  sections  the  facility  recentLy  developed  at  the 
PIBAL,  expressly  for  the  study  of  ablation,  is  described.  A  description 
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of  the  models  and  test  procedure  used  to  obtain  experimentally  the  neces¬ 
sary  information  is  also  included. 


III.  DESIGN  AND  DEVELOPMENT  OF  THE  FACILITY 


The  facility  which  has  been  constructed  for  the  experimental  inves¬ 
tigation  of  ablation  is  depicted  schematically  in  Figs.  9,  10,  and  11.  It  is 
a  supersonic,  blowdown  wind  tunnel  with  an  axisymmetric,  open  jet  test 
section  approximately  six  inches  in  diameter.  This  tunnel  is  connected  to 
the  multiple  outlet  of  the  hypersonic  test  facility  described  in  references  42 
and  43.  A  pebble-bed  heater  with  a  design  temperature  of  3000°R  and  a 
maximum  pressure  of  600  psia  is  the  basic  energy  source  used  in  connection 
with  several  tunnels  and  special  test  rigs. 

The  nozzle  is  a  fixed  geometry,  axisymmetric,  contoured  nozzle, 
which  provides  flow  at  Mach  number  4.4.  Although  fabricated  entirely  from 
type  304  stainless  steel,  the  subsonic  and  throat  regions  of  this  nozzle  are 
cooled  by  a  mixture  of  air  and  water  at  stagnation  pressures  varying  from 
150  to  600  psia,  depending  on  the  nozzle  stagnation  conditions. 

The  contours  for  the  supersonic  portion  of  the  nozzle  were  taken 
from  reference  44.  Corrections  for  boundary  layer  growth  were  made 
according  to  the  method  developed  by  Ruptash  (reference  45)  to  compute  the 
displacement  thickness.  The  nozzle  Mach  number  was  chosen  to  be  4.4  in 
order  that  the  tumelxould be  exhausted  directly  to  the  atmosphere.  Although 
the  PIBAL  facility  has  available  a  vacuum  system  capable  of  reducing  back 
pressure  to  less  than  i/2  psia,  it  was  decided  not  to  exhaust  to  this  system 
during  ablation  tests  because  of  the  possible  generation  of  corrosive  products. 
Hence,  since  the  stagnation  pressure  has  an  upper  limit  (i.e.  ,  600  psia),  the 
Mach  number  chosen  was  that  believed  to  be  the  maximum  for  which  the 
tunnel  could  be  started  with  a  crude  diffuser  and/or  ejector. 

The  arrangement  of  the  test  section  is  shown  schematically  in 
Fig.  10.  An  open  jet  configuration  was  utilized  so  that  the  ablating  models 
could  be  introduced  into  and  retracted  from  the  air  flow  at  well  defined 
times.  The  time  interval  required  for  the  model  to  reach  its  fully  inserted 
position  is  on  the  order  of  1/10  second.  A  pneumatic  cylinder  activated 
the  model  support  system. 
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FIG.  9  GENERAL  SCHEMATIC  -  MACH  4.4  WIND  TUNNEL 


SCHEMATIC  OF  TEST  SECTION 

FIG.  10.  DETAILS  OF  TEST  SECTION  INSTALLATION 


TEST  SECTION  -  CALORIMETER  IN  POSITION 
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FIG.  II  SCHEMATIC  OF  OPTICAL  SYSTEM 


The  optical,  system  shown  in  Fig.  11  utilizes  two  12-inch  diameter 
parabolic  mirrors  with  4-foot  focal  Length  for  schlieren  photographs,  while 
for  shadowgraph  pictures  only  one  is  used  as  shown  in  Fig.  11.  The  test 
section  is  provided  with  10-inch  diameter  windows  of  one-inch  thick,  good 
(but  not  optical)  quality  plate  glass  set  in  steel  frames.  These  are  tilted  at 
an  angLe  of  5°  to  the  axis  of  the  mirrors  to  eliminate  extraneous  reflections. 

It  is  to  be  noted  that  frosted  glass  is  added  to  the  window  nearest  the  camera 
when  the  shadowgraph  technique  is  used.  The  knife  edge,  light  source  and 
cameras  are  installed  outside  the  protective  enclosure.  Adjustment  and 
alignment  of  the  optical  system  can  be  accomplished  externally. 

It  was  originally  assumed  that  sufficient  diffusion  could  be  provided 
by  an  axisymmetric  diffuser  and  considerable  effort  was  expended  to  make  a 
suitable  diffuser.  Starting  was  achieved  for  an  empty  test  section  configura¬ 
tion  with  an  area  ratio  (diffuser  throat  to  test  section  area)  of  1.26.  However, 
tunnel  start  with  a  three-inch  diameter  blunt  model  was  not  achieved.  It  was 
found  necessary  to  use  ejectors  to  reduce  back  pressure. 

The  ejector  and  mixing  tube  are  shown  schematically  in  Fig.  9. 

With  a  throat  area  of  0.88  in.  and  with  supply  conditions  Tstag=540°R, 
Pstag=2000  psia,  approximately  34^/sec.  of  mass  flow  is  provided.  Tests 
have  indicated  that  the  test  section  pressure  is  reduced  by  about  50  percent 
by  this  ejector,  and  starting  with  the  aforementioned  model  has  been  suc¬ 
cessfully  achieved  for  heater  pressures  greater  chan  500  psia.  Detailed  dis¬ 
cussion  of  this  type  of  starting  technique  may  be  found  in  reference  46. 

An  additional  feature  of  the  exhaust  system  is  the  water  spray  for 
cooling  and  purging  of  the  heater  air  prior  to  exhaust  to  the  atmosphere.  The 
ejector  air  also  provides  additional  cooling. 


IV.  DESCRIPTION  OF  MODELS, 

TEST  PROCEDURE  AND  INSTRUMENTATION 


A.  Description  of  Models 

The  parameters  which  are  of  interest  in  ablation  tests  have  been 
described  in  Section  II.  For  the  experimental  determination  of  these  quan¬ 
tities  two  types  of  models  have  been  designed  and  fabricated.  The  first  of 
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these,  which  will  be  referred  to  as  the  calorimeter,  is  shown  in  Figs.  12 
and  14....  This  model  is  used  to  determine  experimentally  the  zero  ablation 
heat  transfer  rates  { c Jq)  and  the  pressure  distribution.  The  calorimeter  was 
fabricated  from  number  type  304  stainless  steel  with  a  wall  thickness  of 
0.050  inches;  it  was  instrumented  with  surface  pressure  taps  and  iron- 
constantan  thermocouples  located  as  shown  in  Fig.  12.  To  insure  that  the 
thermocouples  read  the  surface  temperature  of  the  calorimeter,  the  instal¬ 
lation  was  accomplished  in  the  following  manner:  A  0.040  inch  diameter 
hoLe  was  drilled  through  and  perpendicular  to  the  surface  of  the  model;  the 
thermocouple  wire  was  passed  through  the  hole  and  held  in  place  by  driving 
a  stainless  steel  tapered  pin  into  the  hole  from  the  outside,  crushing  the 
wires  together  at  the  outer  surface  of  the  model.  The  pin  was  then  filed 
flush  with  the  model  surface.  The  external  geometry  of  the  calorimeter  is, 
of  course,  identical  with  the  plastic  ablators  which  are  described  below. 

The  thickness  chosen  is  consistent  with  the  transient  technique  which  is 
utilized  to  determine  the  surface  heat  transfer  rates. 

The  geometry  chosen  for  the  ablating  specimens  was  based  on  the 
following  requirements:  that  the  models 

(a)  be  large  enough  to  permit  measurements  to  be  obtained 
around  an  extended  surface  and  to  provide  transitional 
and  turbulent  flow; 

(b)  have  a  configuration  which  would  be  of  practical  interest 
and  at  the  same  time  be  of  the  simplified  type  for  which 
detailed  analysis  is  available; 

(c)  be  sufficiently  thick  for  a  semi -infinite  slab  analysis  to 
apply,  and  for  the  installation  of  temperature  sensing 
devices  below  the  surface. 

The  final  configuration  chosen  is  shown  in  Figs.  13  and  14  ;  also 
shown  are  the  location  and  depth  of  the  thermocouples  which  are  imbedded 
in  the  plastic.  It  will  also  be  noted  that  a  stainless  steel,  cylindrical, 
afterbody  has  been  added  to  the  model  mount.  This  skirt  has  been  instru¬ 
mented  with  surface  thermocouples  identical  to  those  on  the  calorimeter. 

It  should  be  noted  that  the  values  of  q0  thus  obtained  correspond  to  the 
zero  ablation  case  since  they  are  evaluated  with  the  surface  temperature 
equal  to  the  initial  temperature  at  which  time  no  ablation  of  the  plastic 
models  has  occurred. 

The  installation  of  the  thermocouples  on  the  plastic  models  was 
accomplished  as  follows:  a  series  of  0.125-inch  diameter  holes  were  drilled 
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FIG.  13  DETAILS  OF  ABLATING  MODELS  AND  INSTRUMENTATION 
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through  the  plastic  models  in  the  radial  direction.  Plugs  of  approximately 
0. 125-inch  diameter  of  material  identical  to  that  of  the  model  were  then 
fitted  with  iron-constantan  thermocouples  at  the  desired  depth  as  shown 
in  Fig.  15.  These  plugs  were  thereupon  pushed  into  the  holes  in  the  model 
from  the  rear  or  concave  side  till  the  external  ends  could  be  machined  to 
match  the  hemispherical  surface  of  the  models.  Due  to  the  natural  flex¬ 
ibility  of  the  plastic  materials  the  installation  was  essentially  a  "force  fit" 
(i.  e. ,  the  diameter  of  both  the  plug  and  the  hole  was  essentially  the  same) 
and,  as  a  result,  a  very  smooth  model  surface  was  obtained  (see  Fig.  14). 

Thirteen  plastic  models  were  fabricated,  identical  in  size  and  shape 
and  differing  only  in  the  material  which  was  used.  A  list  of  these  is  pro¬ 
vided  in  Table  I  along  with  data  concerning  the  type  of  material,  numbering, 
etc. 


Some  comments  concerning  requirement  (a)  above  would  be  appro¬ 
priate  here.  Several  investigators  have  found  (see  references  4  and  47) 
that  for  values  of  Nr^>300  transition  will  occur.  This  criterion  was  used 
to  select  the  diameter  of  the  models.  Preliminary  calculations  based  on  a 
hemispherical  geometry  and  a  Newtonian  pressure  distribution  using  Eq.  (43) 
indicated  that  for  a  model  of  3"  diameter  a  value  of  Nr^  equal  to  300  would 
occur  at  a  latitude  of  approximately  60°.  Fig.  lb  shows  the  actual  variation 
of  Nr^  with  "s  as  obtained  from  the  experimentally  determined  pressure 
distribution  on  a  3"  diameter  hemisphere  (Fig.  17). 


B.  Test  Procedure  and  Tunnel  Conditions 

The  experimental  procedure  is  as  follows!  With  all  recording 
potentiometers  running  and  the  heater  stagnation  pressure  and  temperature 
at  the  desired  values,  hypersonic  flow  is  established  by  first  activating  the 
ejector  and  then  withdrawing  the  valve  in  the  throat  of  the  nozzle.  The 
plastic  model  is  then  rotated  into  the  stream;  this  motion  activates  a 
motion  picture  camera  and  an  electric  clock,  so  that  the  model  and  clock 
are  photographed  simultaneously  and  continuously  from  the  instant  of 
exposure.  The  model  is  withdrawn  prior  to  tunnel  shutdown,  the  total  time 
of  exposure  thereby  being  recorded.  The  procedure  is  identical  in  the  case 
of  the  heat  transfer  runs  with  the  calorimeter  model  except  that  motion 
pictures  are  not  necessarily  taken. 

Test  running  times  are  sufficient  to  achieve  steady  state  ablation 
on  a  major  portion  of  the  model.  (Note  that  because  of  the  varying  local 
values  of  qQ  the  steady  state  ablation  rate  occurs  at  varying  times  on  the 
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model.)  These  times  are  estimated  to  be  from  10  to  15  seconds  (see,  for 
example.  Fig.  4).  Table  II  summarizes  the  test  section  conditions  to  which 
the  models  were  exposed. 


C.  Instrumentation 

The  pressures  on  the  surface  of  the  calorimeter  are  measured 
during  the  test  on  transducers  and  recorded  on  single-channel  recording 
potentiometers.  The  pitot  pressure  and  stagnation  temperatures  are  meas¬ 
ured  on  probes  mounted  on  brackets  attached  to  the  end  of  the  nozzle.  A 
schematic  of  the  temperature  probe  as  well  as  a  photograph  are  shown  in 
Figs.  18  and  19.  All  of  the  model  thermocouples  are  also  recorded  on 
recording  potentiometers. 

Because  of  the  open  jet  configuration  which  is  used,  shielding  of 
alL  the  sensing  leads  is  required.  This  is  complicated  by  the  motion  of  the 
model  support  which  is  necessary  for  retraction  of  the  model.  The  shielding 
was  ultimately  accomplished  by  passing  the  leads  through  a  stainless  steel, 
pressure  tight  flexible  hose  as  manufactured  by  the  Titeflex  Company. 

Fig.  20  shows  a  model  mounted  on  its  support  to  which  is  also  attached  the 
aforementioned  flexible  hose. 

The  clock  (see  Fig.  21),  which  provides  time  resolution  of  the  sur¬ 
face  recession,  was  especially  built  for  this  purpose.  For  the  drive  a 
synchronous  motor  with  a  speed  of  2  r.p.m.  was  utilized  in  combination 
with  bevel  gears  with  a  2: 1  gear  ratio  yielding  a  clock  speed  of  4  r.p.m. 

A  back -lighted,  frosted  glass  dial  face  was  provided,  scribed  with  120  divi¬ 
sions.  Thus,  time  intervals  can  be  read  directly  in  eighths  of  a  second. 

The  camera  which  is  used  for  taking  motion  pictures  during  ablation 
runs  is  a  16  mm.  Bell  and  Howell  "Filmo,"  type  70-H.  It  is  operated  by  a 
variable  speed  electric  motor  which  will  give  from  8-64  exposures  per 
second.  A  camera  speed  of  64  frames  per  second  has  been  used  with 
acceptable  results.  The  lens  which  is  used  is  a  Pan-Cinor,  wide  angle  lens 
f2.8.  A  Tiffen  series  8  -  ft  2  Portra-lens  is  utilized  in  conjunction  with  it, 
giving  a  field  diameter  of  10  inches  at  a  subject -to -lens  distance  of  24 
inches. 


The  clock  and  camera  are  activated  from  a  pressure  switch  con¬ 
nected  to  the  activating  cylinder  of  the  model  support,  so  that  the  clock  and 
camera  Start  as  the  model  is  exposed  to  the  flow. 
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V.  DATA  REDUCTION 


A.  Pressure  Distribution 

The  pressure  distribution  which  was  obtained  from  test  runs  with 
the  calorimeter  is  presented  in  nondimensionaL  form  in  Fig.  17.  The  results 
are  in  good  agreement  with  the  modified  Newtonian  theory  of  reference  48. 
Calculation  of  the  theoretical  heat  transfer  distributions  which  are  utilized 
later  are  based  on  these  experimental  values. 


B.  Heat  T ransfer  Rates 

The  local  heat  transfer  rates,  q^,  are  obtained  from  the  surface 
temperature  history  by  the  transient  technique  described  in  reference  49. 

Here  the  subscript  i  indicates  that  the  calculated  values  correspond  to  a 
surface  temperature  equaL  to  its  initiaL  (i.e.,  before  exposure)  constant 
value,  T^.  The  data  is  presented  in  terms  of  the  Nusselt  number  Njqu  and 
is  plotted  against  the  nondimens ional  surface  coordinate's  in  Figp.  22-25.  The 
results  are  compared  to  the  Laminar  and  transitional  heat  transfer  theory 
discussed  in  Section  II,  C.  of  this  report.  In  Fig.  26 there  is  presented 
similar  data  obtained  from  the  cylindrical  after  bodies  during  the  ablation 
runs.  The  calculation  for  the  cylindrical  portion  is  based  on  the  pressure 
distribution  of  reference  56  obtained  for  a  body  with  an  identical  configura¬ 
tion  but  at  a  flow  Mach  number  of  6. 


C.  Ablation  Runs 

Of  the  Thirteen  plastic  models  which  were  fabricated,  eight  (two 
each  of  four  different  plastics)  were  exposed  to  heated  air  for  times  ranging 
from  9  to  28  seconds.  The  tunnel  conditions  and  running  times  for  each 
model  are  presented  in  Table  II  along  with  other  pertinent  data. 

Although  shadowgraph  motion  pictures  were  taken  of  the  abLating 
models,  poor  resolution  of  the  resulting  film  did  not  permit  an  Accurate 
determination  of  the  time  dependent  surface  recession  s(t),  and  thus  also  of 
the  steady  ablation  velocity  B  or  mass  ablation  rate  m  (m=pB).  However, 
the  final  configuration  of  the  ablators  was  measured  at  severaL  azimuthal 
planes  with  considerable  precision  on  a  Jones  and  Lamson  optical  com¬ 
parator  at  a  magnification  of  ten.  It  was  found  that  the  models  ablated  in 
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FIG.  22  DISTRIBUTION  OF  NUSSELT  NUMBER 
ON  CALORIMETER.  RUN  NO.  4  N„«  1.20X10* 


FIG.  23  DISTRIBUTION  OF  NUSSELT  NUMBER 
ON  CALORIMETER.  RUN  NO.  5  NR-I.35XI0* 
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FIG.  24  DISTRIBUTION  OF  NUSSELT  NUMBER 
ON  CALORIMETER.  RUN  NO.  6  H,  - 124  X  I <3* 


FIG.  25  DISTRIBUTION  OF  NUSSELT  NUMBER 
ON  CALORIMETER.  RUN  NO.  7  NR- 1.24X10* 
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an  axisymmetric  manner.  The  resuLts  of  these  measurements  are  pre¬ 
sented  in  Figs  .  27-34  in  absolute  terms  and  in  a  normalized  form  (percent  of 
stagnation  point  ablation). 

The  total  time  of  exposure  t*  was  read  from  both  the  films  and  the 
temperature  traces  and  it  is  assumed  that  the  total  surface  recession  st 
which  was  measured  corresponds  to  this  time  interval;  that  is,  for  t  >  tj- 
no  removal  of  material  occurs. 

The  temperature  histories  obtained  from  the  thermocouples  located 
below  the  surface  are  presented  in  a  semi -logarithmic  plot  in  Figs  35-39.  As 
has  been  noted  in  Section  U,  by  presenting  the  data  in  this  form  the  steady, 
linear  variation  of  the  temperature  history  may  be  observed  from  which  an 
experimental  determination  of  the  quantity  Ba/x  is  obtained. 

We  may  note  at  this  time  that  since  the  steady  mass  ablation  rate 
m  was  not  determined,  it  is  not  possible  to  calculate  the  experimental 
values  of  heff  for  the  various  materials.  Because  of  the  initial  transient 
period,  the  quantity  B^At  *8  ^eBS  than  B,  since  ds/dt^B.  Evidently,  the 
greater  the  period  of  unsteady  ablation,  the  greater  will  be  the  difference 
between  these  two  quantities.  For  a  qualitative  comparison  between  the 
various  materials  a  rough  value  of  heff  has  been  calculated  for  the  stagna¬ 
tion  point  where  the  transient  period  is  a  minimum  in  most  cases  by  taking 
the  measured  ratio  *tAt  as  being  equivalent  to  the  steady  ablation  velocity. 
These  values  are  included  in  Table  II. 

Photographs  of  the  final  configuration  of  the  models  superimposed 
on  a  silhouette  of  the  original  unablated  models  may  be  found  in  Figs.  40-43. 


WADD  TN  60-273 


51 


FIG.  27  DISTRIBUTION 
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FIG.  30  DISTRIBUTION  OF  TOTAL  ABLATION  a*  -  MODEL  L9 
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FIG.  32  DISTRIBUTION  OF  TOTAL  ABLATION  -  MODEL  P6 


FIG.  33  DISTRIBUTION  OF  TOTAL  ABLATION  s,  -  MODEL  N3 


FIG.  34  DISTRIBUTION  OF  TOTAL  ABLATION  st  -MODEL  N7 


FIG  .  35  TEMPERATURE  HISTORY  BELOW  ABLATING  SURFACE  -  MODEL  T5 


o 


<v 

i)x]ui 

*1 


i 

! 


WA DO  TN  60-171 


62 


FIG.  37  TEMPERATURE  HISTORY  BELOW 


FIG.  38  TEMPERATURE  HISTORY  BELOW  ABLATING  SURFACE  -  MODEL  L9 
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FIG.  41.  FINAL  CONFIGURATION  OF  ABLATED  MODELS 


FIG.  42.  FINAL  CONFIGURATION  OF  ABLATED  MODELS 


VI.  COMPARISON  WITH  THEORY  AND  DISCUSSION  OF  RESULTS 


A.  Heat  Transfer  Rates  from  the  CaLorimeter  Models 

Inspection  of  the  data  presented  in  Figs.  22-25  indicates  the  possible 
existence  of  transitional  heat  transfer  on  the  downstream  portions  of  the 
calorimeter.  In  order  to  affect  a  comparison  with  theoretical  predictions 
the  location  of  transition  was  estimated  from  these  plots.  Then,  based 
on  the  pressure  distribution  of  Fig.  17,  the  transitional  theory  described 
previously  was  utilized  to  obtain  the  corresponding  distribution  of  Nj^ju 
The  resulting  theoretical  calculations  of  the  combined  laminar  and  transi¬ 
tional  distribution  are  in  reasonable  agreement  with  the  experimental  data 
as  can  be  seen  in  Figs .  22-25.  It  is  interesting  to  note  that  for  Run  No.  5  with 

Nr  =  1.47xl06  transition  occurs  approximately  at.  a  value  of  "s  =  0.7  whereas 
^  6 

for  Run  No.  4  with  Nr  =  1.20x10  transition  occurs  further  downstream  at  a 
value  of  "s=0.8.  The  corresponding  values  of  Np_^  are  approximately  equal 
and  have  the  value  of  230.  This  is  somewhat  lower  than  the  criterion  for 
transition  which  was  cited  previously. 


B.  Ablation 

For  any  particular  ablation  run  the  tunnel  conditions  and  the  physical 
properties  of  the  model  are  specified,  and  hence  Eqs.  (24)  and  (25)  may  be 
utilized  to  obtain  the  theoretical  variation  of  the  surface  recession  and 
temperature  (at  a  specified  depth)  with  time.  In  particular,  there  may  be 
computed  st=s(tt)  for  several  points  on  the  hemisphere;  these  values  will  of 
course  differ  because  of  the  different  values  of  the  aerodynamic  heat  flux 
qQ  which  exist  at  these  points.  Evidently,  in  order  to  carry  out  the  calcula¬ 
tion  the  distribution  of  qQ  must  be  specified.  It  has  been  already  shown  that 
this  distribution  will  probably  be  transitional  in  character  and  therefore 
highly  sensitive  to  the  particular  flow  Reynolds  number  which  exists.  Since, 
in  general,  the  Reynolds  number  of  any  particular  test  was  not  equal  to  those 
obtained  in  the  calorimeter  runs,  the  following  approach  was  adopted  for 
determining  the  appropriate  heat  transfer  rate  distributions:  From  inspec¬ 
tion  c£  the  data  in  Figs.  27-34,  a  transition  point  was  estimated;  the  heat  transfer 
distribution  was  then  computed  in  a  manner  exactly  analogous  to  that  des¬ 
cribed  previously.  These  distributions  are  included  in  Figs.  27-34  in  a 

Actually,  several  different  values  were  used  until  "best  fit"  of  the  data  was 

obtained. 
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normalized  form  and  are  used  to  carry  out  the  numerical  calculation  of  the 
transier-t  ablation  theory.  Qualitatively  at  least,  we  note  that  the  final 
configuration  of  all  models  corresponds  rather  well  to  the  assumed  (Jq  dis¬ 
tribution;  i.  e.  ,  the  greatest  removal  of  material  occurs  at  stations  where 
the  heating  is  a  maximum.  In  some  cases  of  course,  good  quantitative 
agreement  is  also  in  evidence  as,  for  example,  for  models  T10  and  L8. 

Before  proceeding  to  a  discussion  of  the  results  and  the  comparison 
with  the  theoretical  predictions  it  is  necessary  to  make  some  comments 
concerning  the  properties  of  the  various  plastics  which  were  tested.  In 
particular,  the  materials  which  were  used  fall  into  two  classes;  those 
which  melt  (Nylon  and  Polyethylene),  and  those  which  sublime  (Teflon  and 
Lucite).  Thus  the  calculations  must  be  carried  out  for  this  distinctly  dif¬ 
ferent  behavior  and  the  appropriate  formulation  as  developed  in  Section  II 
must  be  utilized. 

For  the  subliming  materials  Eqs.  (13),  (24),  and  (25)  are  used  to 
determine  s*  and  the  temperature  history  T(x0,t)  where  Xq  represents  the 
depth  at  which  the  thermocouples  were  installed.  Here  the  value  of  II  is 
given  by  Eq.  (1-3)  in  Appendix  I .  The  results  are  presented  in  Figs.  27«30 
and  35-_38.  With  the  exception  of  model  T5,  the  agreement  with  the  experi¬ 
mental  data  is  quite  good,  not  only  at  the  stagnation  point  but  ov^r  the 
entire  hemispherical  surface.  On  the  basis  of  these  results  it  would  appear 
that  the  vapor  injectioxi  at  the  stagnation  region  did  not  influence  to  any 
significant  degree  the  conditions  on  downstream  portions  of  the  model.  The 
rather  large  discrepancies  which  are  observed  for  the  case  of  model  T5 
may  be  attributed  to  the  very  small  amounts  of  materials  which  were 
ablated,  with  the  consequent  increase  in  any  error  due  to  the  measuring 
technique  used. 

Consider  now  the  temperature  histories  which  are  presented  in 
Figs.  35-38.  As  may  be  seen  in  Fig.  2  at  t=tB  the  Eqs.  (6)  and  (13)  give  rise 
to  a  discontinuity  in  the  temperature  history  at  any  given  value  of  x.  This 
behavior  is  of  course  due  to  the  assumed  form  for  the  temperature  profile. 
Thus  it  is  inconsistent  to  compare  the  actual  theoretical  values  of  the 
temperature  at  any  time  with  those  obtained  experimentally.  However,  it 
is  consistent  and  of  interest  to  compare  the  steady  state  slopes  since  the 
integral  method  which  has  been  formulated  does  yield  the  exact  value  of  this 
slope,  as  has  been  indicated  previously.  In  Figs.  35-38  this  comparison  is 
shown  where.for  convenience;  the  theoretical  curve  has  been  shifted  parallel 
to  the  ordinate  (i.e.,  downward) by  an  arbitrary  amount.  Again  fairly  good 
agreement  is  observed,  particularly  for  models  T5  and  T  10. 

For  the  materials  which  vaporize,  reaction  of  the  vapor  with  the 
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constituents  of  the  external  boundary  Layer  may  occur.  It  has  been  indi- 
cated  previously  that  such  combustion  wiLL  have  a  ngeligible  effect  for 
sufficiently  high  values  of  the  driving  enthaLpy  h^-hg.  In  the  teBts  per¬ 
formed  here  this  condition  is  not  fuLfiLLed,  so  that  if  combustion  does  occur, 
the  ablation  rates  shouLd  be  much  greater  than  predicted  by  the  transient 
theory  to  which  the  experimental  resuLts  are  compared.  The  results  of 
this  comparison  would  therefore  indicate  that  combustion  did  not  occur, 
either  for  the  TefLon  or  Lucite  models.  The  physical  evidence  avaiLabLe, 
though  not  conclusive,  appears  to  verify  this  conclusion.  No  combustion  is 
indiated  either  by  the  motion  pictures  of  the  test  runs  or  by  visuaL  inspec¬ 
tion  of  the  ablated  models. 

Consider  now  the  calculations  for  the  meLting  materials.  These 
are  carried  out  as  before,  except  that  now  we  set  n=0  in  Eqs.  (14)  and  (20). 
The  resuLts  of  these  calculations  (see  Figs.  31-34)  yield  values  for  the  total 
surface  recession  s^  which  are  very  much  Less  than  those  actually  exper¬ 
ienced  by  the  ablating  models.  This  result  is  surprising  in  view  of  the 
conservative  nature  of  the  meLting  analysis  which  has  been  utilized  (see 
Section  II).  One  would  expect  this  theory  to  overestimate  the  amount  of 
material  removed,  particularly  at  the  stagnation  point  where  no  build-up  of 
liquid  and  consequent  solidification  can  occur. 

One  possible  explanation  for  this  behavior  is  that  the  values  of  the 
physical  constants  used  for  nylon  and  polyethylene  are  in  error.  In  gen¬ 
eral,  the  problem  of  obtaining  reliable  data  on  the  properties  of  ablating 
materials  is  difficult.  Of  all  of  the  data  which  has  been  obtained  for  use  in 
these  experiments  the  most  doubtful  were  the  Latent  heat  of  fusion  L  and  the 
melting  tempe  rature  Tm  for  the  above-named  materials.1*1  Calculations 
indicate  that  the  value  of  s^  is  indeed  sensitive  to  these  parameters.  (It  is 
also  sensitive  to  (fo  to  the  same  degree,  but  in  view  of  the  consistent  results 
obtained  for  the  subliming  cases  there  is  no  reason  to  suppose  that  the  dis¬ 
crepancy  is  due  to  the  use  of  unrealistic  values  of  the  heat  flux  rates).  By 
adjusting  the  vaLues  of  L  and  Tm  an  attempt  to  fit  the  experimental  data  was 
made.  It  was  found  that  a  unique  set  of  values  for  these  parameters  could 
be  obtained  which,  when  combined  with  the  transient  theory,  gave  reasonable 


We  may  note,  for  example,  that  one  source  (reference  50)  gives  for  the 
melting  temperature  of  nylon  a  value  which  is  higher  than  the  stagnation 
temperature  to  which  the  nylon  specimens  were  exposed  in  this  investi¬ 
gation.  This  would  indicate  that  no  ablation  should  occur.  Nevertheless, 
a  considerable  amount  of  materiaL  was  removed,  and  t-us  removal  was 
evidently  not  a  mechanical  one  as  could  be  seen  by  visuaL  inspection  of  the 
ablated  modeLs. 
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agreement  with  the  experimental  results  for  both  tests  of  each  material. 
These  calculations  are  not  presented  since  it  is  felt  that  the  quantity  of 
data  available  is  insufficient  to  justify  any  definite  conclusions. 

In  view  of  the  possible  inaccurate  knowledge  of  the  physical  param¬ 
eters,  only  one  plot  of  the  temperature  history  is  presented  and  compared 
with  the  theoretically  predicted  values  (see  Fig.  .39).  Again  the  discrep¬ 
ancy  is  evident  and  indicates  the  same  qualitative  trend;  that  is,  the  theory 
gives  a  more  shallow  slope  which  corresponds  to  a  lower  value  of  the 
steady-state  ablation  rate.  This  behavior  is  observed  in  all  cases  for  the 
melting  material* . 


VII.  CONCLUDING  REMARKS 


An  experimental  and  theoretical  investigation  of  the  transient 
ablation  of  plastic  models  in  a  hypersonic  wind  tunnel  has  been  carried  out. 
The  existing  theoretical  studies  of  ablation  have  been  reviewed.  The 
approximate  analysis  of  Goodman  has  been  modified  to  include  the  inter¬ 
action  between  the  vaporizing  material  and  the  external  boundary  layer. 

An  experimental  rig  permitting  tests  in  a  free  jet  at  a  Mach  number 
of  4.4  has  been  designed,  constructed  and  connected  to  the  pebble  bed 
heater  of  the  PIBAL  hypersonic  facility. 

The  models  tested  were  hemispherical  with  a  diameter  of  three 
inches.  The  models  were  instrumented  with  thermocouples  initially  below 
the  surface  and  were  observed  photographically  during  exposure  to  the  air 
stream.  A  calorimeter  model  was  tested  so  as  to  provide  the  heat  transfer 
without  mass  transfer.  Laminar  and  transitional  heat  transfer  were  found 
to  exist  on  the  models;  as  an  interesting  result  of  the  study  it  was  found 
that  the  existing  theories  for  transitional  heat  transfer  under  hypersonic 
conditions  could  be  simply  modified  to  provide  better  agreement  between 
theory  and  these  and  other  experiments. 

The  predicted  total  mass  removed  and  the  steady  state  ablation 
rates  were  found  to  be  in  good  agreement  with  the  measured  data  for  the 
subliming  materials  (Teflon  and  Lucite).  However,  for  the  materials 
which  melted,  the  material  removed  was  found  to  be  considerably  less  than 
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that  predicted  by  the  analysis.  There  are  indications  that  this  discrep¬ 
ancy  is  due  to  a  Lack  of  satisfactory  knowledge  of  the  physico-chemical 
properties  of  the  materials. 
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APPENDIX  I 

DEPENDENCE  OF  HEAT  FLUX  ON  ABLATION  RATE 


Eq.  (14)  indicates  that  the  net  heat  flux  to  the  surface  decreases 
linearly  with  the  ablation  rate.*  This  effect  is  exactly  analogous  to  so-called 
mass  transfer  cooling  resulting  from  vapor  injection  into  the  boundary  layer. 
In  this  case  the  vapor  injection  is  provided  by  the  sublimation  of  the  solid 
material.  Extensive  investigation  of  this  phenomenon  (see,  for  example, 
reference  51)  has  resulted  in  the  following  approximate  correlation  formula: 


R=  9o  - 


<hco'hw) 


(1-1) 


Here  qc  is  the  heat  transfer  rate  to  the  surface  in  the  absence  of 


The  ablation  rate  is  actually  (p  — ^-) ,  but  since  the  density  of  the  plastics  is 

dt 

assumed  constant  in  this  analysis  the  distinction  is  irrelevant  to  the  dis¬ 
cussion  here. 
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mass  transfer  (evaluated  at  the  surface  temperature  Tw),  M  Is  the  ratio  of 
molecular  weights  of  air  to  the  injected  gas,  and  (hoo-hw)  reprejienta  the 
enthalpy  difference  across  the  boundary  layer.  The  numerical  factor  \  is 
valid  only  for  a  laminar  boundary  layer,  and  would  be  significantly  different 
for  the  turbulent  case  (see  references  52,  53,  and  54).  These  references 
indicate  that  the  factor  .2  should  be  replaced  by  .2  . 


If  Eq.  (I- 1)  is  to  be  applied  at  an  ablating  surface,  hw  and  must 
be  evaluated  at  the  ablation  temperature  Ts#  while  P  and  M  must  correspond 
to  the  gas  liberated  by  the  ablation  process.  By  virtue  of  the  above  con¬ 
siderations  Eq.  (14)  becomes 


so  that 


i 

q(t)  =  qol  -  |  M4p(hQO-h8)  # 

T  =T 

1  x  m 


n 


M,P(hoo-hs> 


(1-2) 


(1-3) 


It  shouLd  be  noted  here  that  chemical  reactions  and  re-radiation 
from  the  ablating  surface  are  implicitly  excLuded  by  the  expression  for  the 
heat  balance  at  the  interface  [i.  e.,  Eq.  (9)].  For  pure  pLastic  materials 
(and  in  general  for  subliming  types  of  materials),  the  neglect  of  radiation 
is  justifiable,  since  the  maximum  surface  temperatures  wiLL  be  relatively 
Low.  No  such  a  priori  justification  exists  for  the  neglect  of  chemical 
reactions.  However,  at  least  one  investigator  (reference  55)  has  found 
experimentally  that  no  combustion  takes  place  under  hypersonic  conditions. 
The  experimental  results  obtained  here  also  indicate  the  absence  of  any 
chemical  reaction  (see  Section  VI). 
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APPENDIX  II 


TEMPERATURE  PROFILE 


The  temperature  profile  assumed  in  the  previous  analysis,  in  a 
slightly  rearranged  form,  is 


T(|'.yT-t  =  exp[-Sjt] 
s  i 


(H-l) 


8  =  s(t) 

6  =  6(t) 

The  quantity  8  has  been  referred  to  at  times  as  the  "e-folding"  length  (see 
reference  21),  since  at  the  distance  x=  8+  s  from  the  origin  the  temperature 
is  given  by 


KS+s.tJ-T  1 


T  -T- 

s  1 


<H -2) 


Of  significance  from  the  mathematical  point  of  view  is  the  fact  that 

(i)  The  profile  of  Eq.  (II— 1 )  approaches  the  exact  profile  for 
steady  ablation  for  large  time. 

(ii)  The  profile  does  not  satisfy  exactly  the  initial  conditions 
which  are  carried  over  from  the  "heat  sink"  portion  of  the 
problem  [i.e.,  Eq.  (12)J. 

With  regard  to  (ii)  above,  the  following  may  be  noted.  First,  the  profiles 
of  Eq.  (H-l)  and  of  Eq.  (6)  give  for  t=ts  identical  temperature  gradients  at 
the  surface,  i.e.: 


8T(0,t„)  qo 
8x  =  '  k 
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Second,  the  heat  stored  in  the  solid  at  this  time  as  calculated  by  the  two 
profiles  differ  roughly  only  by  25%  [the  exact  profile  (6)  gives  a  lower 
value  of  stored  energy]  for  the  types  of  materials  and  the  environment  to 
which  this  analysis  is  applied.  For  example, 


Q 


stored  (exact) 


fCO 

/  Pc  (T-T.)dx 
0 


=  eyToo-Ti)2Nu 


k 

h 


erfc  X  -  e 


Nu(Nu+2X) 


erfc(Nu+X)dX 


=  PyT^-T.JZNu 


k 

h 


l  l  Nu 

[— ~r - (1-e  erfc  Nu)] 

n/tT  2Nu 


where,  of  course,  Nu  is  evaluated  at  t=tg. 

On  the  other  hand,  for  the  approximate  profile 


stored  (approx 


.)  =  /€°PVT.’Ti,eXp[-6^ldX 


(T  -T. )  k 
=  Pc  'a  1 1  * 

P  lToo-T.)  h 


The  ratio  is  thus 


Qst .  ex.  (Too  ~Tj)  (Tqq  "Tg  ) 


st.  app. 


<VTi)3 


2N  u  Nu* 

-  (l-e™  erfc  Nu)} 


(II -3) 


For  Teflon,  based  on  the  values  of  its  properties  listed  in  Table  I  and  the 
conditions  under  which  the  experiments  were  performed,  the  ratio  given 
by  Eq.  (II-3)  is  approximately  0.75.  Since  the  toal  heat  stored  at  this 
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time  is  in  any  case  small  compared  with  the  total  heat  absorbed  during  the 
entire  ablation  process,  this  error  is  quite  negligible  (see  Fig.  2.  for  com¬ 
parison  of  the  two  profiles). 


WADD  TN  60-273 


82 


^  — 


UNCLASSIFIED  I  UNCLASSIFIED 


AM 


^ - ^ - 


& 

E 


C'N-^ 


0 

M 

E 

M 

§ 


« 

B 

fc 

W 

1 


6  5 

s 

i  o  o 


_  42 


o  c 

•H  *H 

"8  c 

t  © 

a 

'  1 

O  *  £  2  H 
Cm.  ©  42  C  © 

C 1  -S  3  8 

lliSI! 

©  © 

£  9 

1 

©  r-*  L 

©  C.  42  © 

©  X  o 

5  o 

1 

X  O  *Q  ©  © 
■H  C  X  Cm* 
42  X  ©  42  L 

1^*1 
6  in  rl  O 

It 

©  O 

1 

13  S3f 

L,  42  ©  42 

S>  ©  x  n 

j)!  .H 1  £ 

'O  n  >>  o  _ 

(0  ©  ©  r-f  ©  *0 

D  >  xh/:  5 

•H  ©  (?  H  O 


•H  J-*  •H 


£  ©  © 


t  R  8.|jh 

■°1  !fs 


E-a  9  8 

©  ©  3  10  O 

ill! 5?  I 


•tr 


- <£ 


H 

o 

25 

22 


42  tj 

©  o 

r  -f  t« 

S3 


u. 

s\f| 

3  * 

o 

# 


■g .. 

3  *H 

5  5 


8 

4J 

c 

.3 

8  8 

g<S 


O 
O  -*H 

•sa 

>  > 


5 

CN 

ft? 

V  VO 

O  H 

i  ©  VO 
I  C.  V— ' 

V  c~v 
C  <^\ 

<S 


fl 

E 


0 

£ 


o 

z 

22 


M 

W 

M 

fc 

U) 

CO 

a 

o 


3K 

£2 

iiM 


t 

O 

a 

5 

2 


.JE. 


J3  © 
©  T3 


&t&Sf>CQ 

h  o  &  z  ^  q 

^  atnp  w  ^ 

8  J* 

«.  top?  ©  •• 

CQ  Cm  z^H\n 

cn  M  x>  o 

<M  •  £-4  V"  *  rt  VO  c 

o  S>»  CO  ©42  OS 

t  W  O  O  •  »A 

3£^g  g-g^A 

3  ©  25  Q  C  ©  ,  v£) 

l<  O  «J  o  hJH 

•h  o  m  cr  o  *<vo 

V.OHWW  w  pc  ^ 
©  (fl  *£  a,  W  H  c^\ 
O  *-«  Cl,  r'v 

m  m  b:  c  q  w 

©  -h  q 

U  (J  ■<  +>  •  < 

■H-fjt  C  rH  <NJ 

C  E  O  Z  rf  (j  v\  +> 

X  ©  h  42  C  VO  C 

g  co  a c 

*0  *-3  W  0  •  -♦J 

opwo  j*:  c 
..two  a  ©  o 
o  ©Wo  >»  CM  ©  o 
a.  *<  ci  >;  pco 


;:£§  las 

5  ?l  ,*  **•  £ 

Hi 


♦>  fH 

<t  © 

A%t 

tz"  t 

©  «H  ©  O  _ _ 

•HrHXOfcCDTJ©* 

S3  P,  So.  S£  9 

©  ©  ©  ©  o 

u  •  X)  tj  x  x 

42  ©  x  ©  ©  42  O 

14  ©  ©  Ck,  *0  P  -H  -H 
C«  O  rH  tT  ©  ft  b  4- 
P£  EH  * 

™I t 


© 

42 

o  c 

Vh  -H 


dj . . 
81*42 

•  r  4>  <h  C  ©  0.  © 

©  42  C  ©  -H  X  r< 

clU3  S3  *!  . 

m  _ i  i.  m  v 


\i  $”§  J?  t 

1  *g  <n  >>  o 

l|ls3ll 

t  -H  ttf  42  bO 

c 


©  X 
X  o 


1  «v  »■» 

I  M  ©  ►.  O  42 

'  3  S  8  8  «  fc 


o  c 


o 

I 


©  -H  a  O  ri  ©  42 

421kU  «  ©  >  C  C 

CUO©  O  ©  O  © 

Sffl  Cm.  X  CT»  -H  x  -H  • 

O  42  J-  -H 
w  ©  ©  je  i 


I  8£ 

g4  *.  C  §  S  l  iA 

«  r-c  O  X.  W  JO 

^3 


C  ©  X  3 
I  ffl  O  -H  '*4 
rH  -H  t.  O 

|  °-g  SL© 

>  ©  rfl  ©  X 
)  ©  v>  ►  3 


Lt. 

•rf  1.-H 

"  6  8  6  5  |  • 

xoxt'bCw.. 

**  **  6  ® 

2  fe  i.  £  8.1  Ri 

■"?n* 

;  "2  o>  > 

§  8  £ 


^8 


•  o  »  • 

x  i  je  5 


13  IS 


0  8. 


©  *1  I.  *H 

4>  L  ©  C.  © 
»  p.t<  rt  f  2  c 

55  8 


I 

4t- 


...  ..gjjSs 

1 8  ?s 

Ife3^3  Sa  •  „ 

rl  c.  (J  <H  O  ©  rH  Q. 

r  E  -HjjX  fi  fl  J 

5 


’C  J  3  « 

a^cfsl 

©do 

3  2  3  o 

*>  «  -,  ..  . 

3-D  .1|« 
i  S  3  a  2  53  p 

8 3  B5  S 

5  3  ®  5  ®  "S  ^ 

•H  42  c  j5  TJ  jj 


«  42  tj  , 

•  H  P>  t 

©  ©  ©  Q  O 


.  JS 
r?  9  „ 

c.  o  © 


i.  42 


*S 

3  8  2  _ 

U  4>  © 

©  ©  MO©  W 

♦>  X  'll  Q»  ^  CO  ^  ©"Q  ** 

©Hcr*oi4©ci© 


